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SIMULATION OF A TORBOFAN ENGINE FOR EVALUATION 


OF MULTIVARIABLE OPTIMAL CONTROL CONCEPTS 

Kurt Seldner 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 


ABSTRACT 

The development of control systems for jet 
engines requires a real-time computer simula- 
tion. The simulation provides cin effective tool 
for evaluating control concepts and problem areas 
prior to actual engine testing. 

This paper describes the development and 
use of a real-time simulation of the Pratt & 
Whitney FlOC-PW-100 turbofan engine. The simu- 
lation is being used in a multi-variable optimal 
controls research program using linear quadratic 
regulator theory. The simulation is used to 
generate linear engine models at selected oper- 
ating points and evaluate the control algorithm. 

To reduce the complexity of the design, it 
Is desirable to reduce the order of the linear 
model. This paper discusses a technique to re- 
duce the order of the model and compares selec- 
ted results between high and low order models. 

The LQR control algorithms will be pro- 
grammed on an SEL 810B digital computer. This 
computer will control the engine simulation over 
the desired flight envelope. 

INTRODUCTION 

The development and evaluation of controls 
for aircraft propulsion systems requires a dy- 
namic model of the engine and controls . A rep- 
resentative simulation provides the capability 
to predict engine performance over the specified 
flight envelope. The simulation also provides 
an effective tool to evaluate control concepts 
prior to hardware implementation on the actual 
engine. In addition, the simulation can serve 
as an aid in studying problems during the devel- 
opment and testing phases of the program. 

T1 paper discusses methods of simulating 
aircratt propulsion systems with particular 
emphasis on a hybrid simulation of the Pratt & 
Whitney FIDO afterburning turbofan engine. Lewis 
Research Center is Involved in a joint NASA-Air 
Force controls program to evaluate a multivari- 
able control for the FIDO engine. The control 
will be based on linear quadratic regulator (LQR) 
theory. This theory requires that llneai' state 
variable models at various nominal operating 
points be known. This paper discusses a method 
to generate these models from the non-linear 
hybrid FIDO simulation. 


The F-lOO engine mudel contains states that 
may not be measurable or the dynamics may not be 
significant for control design, linear state 
regulator theory requires that all states are 
either measurable or observable. To reduce the 
complexity of the LQR multi-variable controller, 
the order of the linear models was reduced to- 
include several selected states based on engine 
controls experience . A discussion of the methods 
employed to generate and reduce full state linear 
models is presented. 

COMPUTER MODELING OF ENGINES 

Modern simulation techniques are mainly 
limited to digital and hybrid computers. The 
digital computer, due to its accuracy and re- 
peatability, is generally preferred. However, 
real-time analysis cannot be easily accomplished 
with the digital computer. For this reason, the 
hybrid computer which incorporates the advantages 
of both analog and digital computers is ideally 
adapted to real-time simulation of aircraft 
engines . 

A major problem encountered in modeling air- 
craft engines is to develop the bi-variate func- 
tions required to describe the performance of the 
rotating components of the engine. The digital 
computer can be programmed to generate these per- 
formance maps, whereas the analog computer is 
best utilized for calculations, such as integra- 
tion, multiplication and single variable function 
generation. 

Whenever combining digital and analog com- 
puters, for real-time simulations, one must con- 
sider the digital computer up-date time (time 
required to perform the map search and digital 
calculations) that cannot be exceeded if the 
closed loop accuracy and stability are to be main* 
tained over the desired frequency range. 

ENGINE SIMULATION 

To support the multi-variable optimal con- 
trols program, a real time hybrid simulation of 
the Pratt & Whitney FIDO Series I afterburning 
turbofan engine was developed (Ref 1) . Figure 1 
Is a schematic diagram, showing the various com- 
ponents of the engine. The approach used to sim- 
ulate this type of engine is to model each com- 
ponent using basic conservation equations . 



A dl&ita!. ctmulatlon, developed by the 
engine •janufacturt.r (F&W), was used to obtain 
both steady state and transient performance 
characteristics. The bivariate functions, de> 
fining the steady state performance of the com- 
pressor, fan and tubine were obtained from this 
sinilation. These bivariate maps define the 
flow, pressure and enthalpy relationships as 
functions of fan and compressor rotational speeds. 

The functions are stored in the digital 
part of the hybrid computer. The computer is 
capable of handling up to twelve functions of 
six pairs of independent variables. Whenever 
functions are highly non-linear (e.g. compressor 
■aap) , special radial interpolation techniques 
are required. A digital subroutine has been 
developed for this type of function generation 
(Ref 2) . Multi-variable functions of a single 
pair of independent variables can also be 
generated . 

For modeling the FlOO engine, volumes were 
assumed for components either where the gas 
dynamics are important or where algebraic loops 
must be avoided. The energy, momentum and con- 
tinuity equations are solved £(ac each volume. 

The rotor dynamics are computed from the con- 
servation of angular momentum. 

To satisfy the real time requirement , 
several model simplifications and assumptions 
were considered. To compute the temperature 
ratio across the compressor and fan, the 
efficiencies have to be included in the bi- 
variate functions. To eliminate the additional 
map search, the temperature ratios were gener- 
ated as piecewise linear functions of pressure 
ratio. Secondly, the gas properties were 
assumed constant, except for the torque balance 
relations, where the values for specific heat 
were scheduled as functions of compressor and 
fan discharge temperatures. This correction was 
necessary to improve the torque balance along 
the normal operating line. 

MODEL VERIFICATION 

To verify the engine model accuracy over 
the specified flight envelope, the required 
values for the control inputs must be supplied 
to the hybrid simulation. The control inputs 
for the FlOO engine over the non-afterbumlng 
region are main burner fuel flow, exhaust nozzle 
area, inlet guide vane and stator vane position. 
The magnitudes of the control inputs and the 
engine output variables for various flight and 
operating conditions were obtained from the P&W 
digital simulation. This evaluation was per- 
formed to insure that the model was valid. 

Figures 2-S present the steady state hybrid 
results with the values obtained from the digital 
simulation. As can be observed, the comparison 
between the two simulations at the sea level 
static condition is satisfactory. Slight dis- 
crepancies can be observed for fan airflow at 
high altitude condition (figure 3) . This error 
can be attributed to Reynolds number effects 
which were not included in the FlOO Series 1 
simulation. An updated simulation includes this 
effect. 


The dynamic accuracy of the engine simula- 
tion mist be acceptable over the specified fll^^t 
envelope. The four control inputs, as previously , 
defined for steady state operation, were scheduled 
as a function of time between idle and inter- 
mediate power conditions. The idle power condi- 
tion is defined as the minimum engine thrust. * t 
The intermediate pmier condition is the maximum 
engine thrust before entering the engine after- 
burning region. Again, the scheduled control in- 
puts were obtained frutn exercising the FSM 
simulation through transients. The effort was 
restricted to evaluating only the dynamic open- 
loop performance of the hybrid engine simulation 
by scheduling the control inputs instead of 
utilizing the P&W bill of materials controls. 

Figure 6 illustrates the transient responses 
for fan speed and nozzle inlet pressure by ad- 
vancing the power lever from idle to intermediate 
power condition. The results for the hybrid and 
digital models show that the hybrid simulation 
adequately represents the dynamics of the engine. 

To verify the dynamic accuracy of the hybrid 
model, the transient responses of all significant 
output variables were determined over the com- 
plete flight envelope. 

Because the engine simulation must also be 
valid for the afterburning region, both steady 
state and dynamic operation was evaluated at this 
condition. The additional control input, after- 
burner fuel flow, was scheduled as an input to 
the hybrid simulation. To study the transient 
performance, the power lever was advanced beyond 
the intermediate power condition to the full 
afterburning mode. 

The approach used for recording and analyzing 
the transient data and, subsequently, the linear 
models and reduction of these linear models is 
shewn by figure 7. The process of recording, 
reducing and plotting the data was completely 
automated. The link between the hybrid and 
IBM 360 computers is a SEL 810B digital computer, 
which was utilized to record the transient data 
and generate the information on paper tape. A 
360 computer program was utilized to crossplot 
the hybrid transient data with the F&W transient 
data. 

LINEAR MODELS 

After the non-linear engine model evaluation 
was completed, the FlOO hybrid simulation was 
utilized to generate linear models. The linear 
models arc useful in analyzing the performance of 
the engine about an operating point and also for 
the linear quadratic regulator design. 

Full state linear models were generated by 
adding a small disturbance to the nominal value 
of the state and cc .itrol variables . The re- 
sultant deviations «ere recorded by the SEL 610B 
digital computer and then processed on the IBM 360 
to compute the various system matrices (A, B, C, D). 
The procedure used to derive and reduce the state 
variable data is also illustrated by figure 7 . 

Step responses for various control inputs can be 
generated directly on the IBM 360 computer. 
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The results denonstrated by the hybrid and 
linear models for a main burner fuel flow step 
disturbance at the design condition are llius- 
trated by figures 8 > 9 . Some error can te ob- 
served for the final steady state values for ; 
both fan speed and exhaust noszle pressure. 

These differences can be attributed to the non- 
linearities of the engine simulation. An analy- 
sis can be performed to determine the optimum 
values for state and control disturbances to 
restrict operation within the linear range. How- 
ever , about an operating point extremely small 
disturbances will result in objectionable noise 
and recording errors in collecting the linear 
model data. 



( 3 ) 


( 4 ) 


— A22» B21 C^« C2 are sub- 

matrices resulting from the partitioning 
the state vector x 


MODEL REDUCTION 

The primary objective of the FIDO hybrid 
simulation is to test and evaluate the optimal 
multi-variable control algorithm. The linear 
quadratic regulator (LQR) design requires 
knowledge of the system matrices (A, B. C, D) at 
various nominal operating points ^out the flight 
envelope. The linear models, described in the 
Model Verification Section, contain all states 
as defined by the hybrid simulation. Because 
LQR theory utilizes full state feedback (either 
measurable or estimated) , the feedback control 
would be too impractical to develop and implement 
with the actual full state model engine repre- 
sentation. Therefore, it would be desirable to 
reduce the higher order linear model to an 
equivalent lower order model. 


= mxl sub-vcctor 
X2 •* (n-m)xl sub-vector 

Aj^i “ nwm matrix 
Aj ^2 " mx(n-m) matrix 
A21 “ (n-m)xm matrix 
A22 “ (n-m)x(n-m) matrix 

- mxp matrix 
B2 ® (n-m)xp matrix 
Cj = rxm matrix 
C2 = rx(n-m) matrix 

Assume that ^2 Solving eq. 3 for X2 and 

substituting fur X2 into the equation yields 


Axp 

+ Bu 

( 5 ) 

Cxr 

+ Du 

(6) 


The difficulty encountered in the model 
reduction was the selection of appropriate states 
for the lower order model. The final choice was 
based strictly on past experience In selecting 
the significant states for this type of engine. 
The method applied to the FIDO full state linear 
model was to partition the state vector into a 
set of slowly decaying states and a set of 
rapidly decaying states. The dynamics of the 
fast response states were neglected by setting 
the derivatives for these states to zero and 
approximating these states by their steady state 
values. The method is described in detail by 
reference 3. 

The mathematical analysis is as follows: 
Consider the linear system 


X “ 

/u 

+ Bu 

(1) 

y * 

Cx 

+ Du 

(2) 

where 

X 

= state vector nxl 



u 

= control vector pxl 



A = state matrix nxn 
B B control matrix nxp 
y "> output vector rxl 
C ■ output matrix rxn 
D B control matrix rxp 

For the FIDO model reduction, a sub-vector 
consists of states that are to be included in 
the reduced model, whereas a sub-vector X2 in- 
cludes the remaining states. 


where x^ is the reduced system state vector 
and the reduced system matrices are: 

-1 

A = (Aj^^ - Aj^ 2A22^2 i) matrix 

B ' (Bi - A12A22B2) ("ixp) matrix 
£ ■ (^1 ■ C2A22A2I) (rxm) matrix 
D = (D - C2A2JB2) (rxp) matrix ■ 

Equations S and 6 represent the reduced order 
linear model. The system matrices (A, B, C, D) 
c^ ^ n'pi^aced by the reduced system matrices 
(A, B, C, D) . These latter system matrices con- 
tain sufficient information to reconstruct the 
eliminated states and the output vector y. 

The optimum reduction could be achieved by 
retaining only the slowly decaying states and 
eliminating the rapidly decaying states. If a 
reasonable separation exists between low fre- 
quency and high frequency states, then the lower 
order model will be a better approximation to the 
higher order model. The eigenvalues computed from 
the higher order model, derived from the hybrid 
simulation, are not widely spaced so that the 
reduced model is only approximate. In addition, 
a mode analysis indicates interaction between 
states. This implies that eigenvalues cannot be 
easily associated with their respective states. 
Depending on the sufficiency of the reduced model 
approximation, the eigenvalues will differ from 
those of the full state model. 

Figures 10-12 compare selected output vari- 
able step responses for two reduced order with 
full state models. The fan speed was a retained 
state for the reduced order model, whereas the 
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hl^ turbine Inlet temperature was reconstruete^ . 
nmiat la an output quantity, t^ch was rectm* 
structed from the reduced order model, the 
results Illustrate that the seventh order system 
is adequate and represents the full order model 
reasonably well, the' difference between the 
various step responses can be attributed to the 
lack of high frequency Information, this fact 
can be readily observed for both the fan and 
turbine Inlet temperature where the response 
time Is decreased as the model order Is redu«^d. 
For the thrust response, due to the feed^forward 
term Du, an Initial thrust Increase was observed 
for the third order reduced model, the effect 
Is not as noticeable for the full state and 
seventh order models. Sll^tly slower response 
can be observed for the third order model. It 
should be emphasized that the selection of the 
retained states was based on intuition and engine 
control experience. Additional analysis Is be- 
ing performed to establish the final order of 
the reduced model and selection of states. 

From these results it can be concluded that 
the full state model can be adequately repre- 
sented by a third order model . The final deci- 
sion for selecting the retained states depends 
upon the monitoring, protection and control 
problems . 

CONCUJSIONS 

This paper has briefly discussed real-time 
simulation of an augmented turbofan engine. The 
results indicate that the hybrid computer Is an 
effective tool to model engines efficiently and 
accurately. The steady state and transient 
performance of the model compares favorably with 
results obtained from the BSrW digital simulation. 

The model can be utilized in evaluating 
modem control concepts and as a tool in solving 
problem areas during the development and testing 
phases of the pi-ogram. 
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Figure 3. > Comparison of hybrid and digital steady-state 
data for fon airflow at sea level static. 
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Figure 4. - Comparison of hybrid and digital steady-state 
data for compressor speed at altitude * 22 SS5 meters, 
(Mach number "2.5 
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Figure 5. - Comparison of hybrid and digital steady-sta 
data for tin airflow at altitude ■ 22 5K meters. Mach 
number -2.5. 
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Rgure 6. - Comparison of hybrid and digital simulation 
transient performance for Idle- to-intermediate power 
lever advance, sea level static conditions. 
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Figure 7. - Analysis of transient data and linear models. 
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Figure 12. - Step responses ^ various linear models for 
fuel flow disturtance. sea level static condition. 


